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i^T-i' ABSTRACT 

o 

Q^) ■ We present the analysis of the large population of star clusters in the blue com- 

pact galaxy (BCG) Mrk930. The study has been conducted by means of a photometric 
analysis of multiband data obtained with the Hubble Space Telescope (HST). We have 
reconstructed the spectral energy distributions of the star clusters and estimated age, 
O . mass, and extinction for a representative sample. Similar to previous studies of star 

clusters in BCGs, we observe a very young cluster population with 70 % of the systems 
formed less than 10 Myr ago. In Mrk930 the peak in the star cluster age distribution 
at 4 Myr is corroborated by the presence of Wolf-Rayet spectral features, and by the 
observed optical and IR lines ratios [O iii]/H/3 and [Neiii]/[Ncii]. The recovered ex- 
tinction in these very young clusters shows large variations, with a decrease at older 
ages. It is likely that our analysis is limited to the optically brightest objects (i.e. 
systems only partially embedded in their natal cocoons; the deeply embedded clusters 
being undetected) . We map the extinction across the galaxy using low-resolution spec- 
tra and the Ha/H/3 ratio, as obtained from ground-based narrow band imaging. These 
results are compared with the extinction distribution recovered from the clusters. We 
CO . find that the mean optical extinction derived in the starburst regions is close to the av- 

■ eraged value observed in the clusters (more than 80 % of systems have E(B — V) ^ 0.2 

mag), but locally, do not trace the more extinguished clusters. Previous HST studies 
of BCGs have revealed a population of young and extremely red super star clusters. 
We detect a considerable fraction of clusters affected by a red excess also in Mrk930. 
The nature of the red excess, which turns up at near-IR wavelengths (/ band and 
longward) remains unknown. We compare the cluster formation and the star forma- 
tion history, the latter derived from the fit of spectral population synthesis models to 
the spectra. We find a general agreement between the two independently estimated 
quantities. Using the cluster properties we perform a study of the host environmental 
properties. We find that the cluster formation efficiency (the fraction of star formation 
happening in clusters) is significantly higher, suggesting a key role of the environment 
for the formation of these massive objects. 

Key words: galaxies: starburst - galaxies: star clusters: general - galaxies: irregular 
- galaxies:star formation 



* Based on observations made with the NASA/ESA Hubble 
Space Telescope, obtained at the Space Telescope Science Insti- 
tute, which is operated by the Association of Universities for Re- 
search in Astronomy, Inc., under NASA contract NAS 5-26555. 
These observations are associated with program # GO 10902. 



1 INTRODUCTION 

Galaxy encounters give rise to vigorous bursts of star 
formation in which numerous of young star clusters are 
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form in g (e.g. Whitmore fc Schweizerl 19951. Holtzman et alj 



19961; lAlonso-Herrero et alj |2000| ; IVaisanen et all l2odsf 
Pellerin et alj|20ld ). However, young star clusters are also 
observed in quiescent spiral galaxies, where the star for- 



mation rate (SFR) is lower and continuos (e.g. Barth et al.l 
19951: ICarlson et al.lll998l ; lLarsenll2002l iBastian et al.ll2005l : 
Mora et al.l l2009h . Massive star clusters form, during 
sporadic episodes, e ven in d warf galaxies (NGC 1569, 
lArp fc Sandagelll985l; N GC1705. iMelnick et ai1ll985l; He2- 
Conti fc Vaccall 19941; SBS 0335-052E. iThuan et al.ll 19971 ; 



10. 



ESQ338 IG-04. lOstlin et~ai]|l99Sl ; iBillett et al.ll2002T i 



The luminous (M s < -18.0), massive (~ 10 9 " 10 M o ) 
blue compact galaxies (BCGs) show clear signatu res of in- 
teractions and/or mergers (e.g. lOstlin et af]|200ll ) and the 
numerous observe d massive clusters are likely formed in 
these encounters (lOstlin et all 120031 ; I Adamo et"ai] l2010al ; 
lAdamo et al]|201 j ). The very bright ultraviolet and optical 
luminosities of these systems suggest rather low dust con- 
tent and metallicity. Spectra dominated by emission lines 
clearly demonstrate that BCGs are undergoing a burst of 
star formation event. The youth of the burst episode is also 
observed in the recovered age distribution of the star clus- 
ters, which shows a peak of cluster formation younger than 
5 Myr. BCGs are considered to be analogs of high redshift 
Lyman break galaxies and can probe galaxy formation and 
evolution at higher redshifts, when the Universe was younger 
and less enriched in metals. 

The analysis of the young cluster populations in BCGs 
is quite challenging due to the rapid evolution a cluster ex- 
periences during the first 10 Myr (still partly in an embed- 
ded phase). Moreover, this analysis is based on the inte- 
grated luminosities of the clusters, which are mostly unre- 
solved at the distance of the targets. Observations of re- 
solved newly born star clusters in the Milky Way and in the 
Magellanic Clouds reveal that these are quite complex sys- 
tems. A cluster forms from the collapse and fragmentation 
of giant molecular clouds (GMCs; Lada & Lada 2003). The 
compactness of the proto cluster determines whether the 
congl omerate of stars will form a cluster or a loose associa- 
tion |Portegies Zwart et alj|201ol ; iGieles fc Portegies Zwartl 
2010). The massive and short-lived stars rapidly reach the 
main sequence and produce strong winds and UV radiation, 
which ionize the intracluster gas and create bubbles and 
shells. These Hn regions surround the optically bright core 
of stars and significantly contribute to the integrated fluxes. 
However, a large fraction of the stars is still accreting mate- 
rial from their dusty disks (young stellar objects, YSOs) or 
contracting (in the pre-main sequence; PMS phase). Sim- 
ple stellar population models assume that PMS stars do 
not contribute to the integrated flux of the cluster. This 
ass umption is valid a t blu er spectral range but not in the 
IR. iMever fc Greissll |2005h estimated that PMS stars con- 
tribute between 5 and 17 % to the total flux in the H band 
during the fi r st 3 M yr of cluster evolution. As follow up, 
iGreissl et all (|2010l ) found direct evidence of PMS stars in 
the spectrum of an unresolved young star cluster in the An- 
tennae syste m. Moreover, the e dges of the clusters are places 
for trigge red (lElmegreen 19981) and progressive st ar forma- 
tion (e.g. lCarlson et al.ll2007l : IWalborn et alj|2002l ). Delayed 
or triggered star formation processes in dense and dusty re- 
gions surrounding the cluster could explain a large fraction 



of massive YSOs contributing to the IR spectrum of a cluster 
a few Myr old. 

A significant fracti on of young star clu sters in Haro 
11 (lAdamo et alj|2010ah and ESO 185-IG13 (|Adamo et al.1 
l201lT shows a clear signatu re of a flux excess in the 
near-IR. The models we use (|Zackrisson et all l200ll . here- 
after Z01) include a self-consistent treatment of the pho- 
toionized gas, important during the first Myr of the 



cluster evolution (e.g. Kriiger et al.l 1995 




Zackrisson et al. 


2001 


; lAnders & Fritze-Alvenslebenl 12003 


Zackrisson et al. 


2008 


; Reines et alj|20ld: Adamo et alj|2010bh. Adamo et al. 



the SEDs of the clusters during the first 10 — 15 Myr of 
cluster evolution. In metal-poor environments, the contri- 
bution becomes smaller after 6 Myr on the blue side of the 
cluster spectrum, but lasts longer in the NIR wavebands, 
contributing between 10-40% o f the total NIR fluxes at ^10 
Myr. Since the models used in I Adamo et"aj] |2010al . l201ll ) 
already include a contribution from photoionized gas, the 
cause of the excess in BCGs star clusters resides in other 
mechanisms (e.g. YSOs, PMS stars, hot dust, etc.). 

Mrk930 is a BCG located at roughly 72 MpcQ. The 
galaxy was imaged and cataloged in the mid-60s and the 70s, 
during the first large-scale objective-prism survey for galax- 
ies wit h blue and ultraviolet excess in their continuum radi- 
ation (|Mazzarella fc Balzanolll986l ). i.e., the First Byurakan 
Spectral Sky Survey. Due to the low meta llicity content 
(12 + log(0/H) = 8.06. llzotov fc Thuadfl99ll . Mrk930 was 
included in the late 90s among the galaxies studied for its 
similarity with h i gh red shift Lyman break galaxies (LBGs). 
Ilzotov fc Thuanl |l998T ) classified Mkr 930 as a Wolf-Rayet 
(WR) galaxy. Typically, WR signatures in the spectrum of 
a low metallicity galaxy is a rare phenomenon because the 
duration of th e WR phase is shor ter when the metallicity 
decreases (see iGuseva et al.1 12000| and references therein) . 
The presence of WR features indicates that the galaxy is 
undergoing a young burst episode. 

Mrk930 was included in Malkan' s Hubble Space Tele - 
scope survey of local active galaxies (jMalkan et al.|[l998l ). 
as a Hn galaxy. The short F606W exposure of the galaxy 
revealed that the two starburst knots observed in ground- 
based images were in rea lity formed by numerous bright 
star clusters |6stlinll2OO0l ). The presence of a very active 
starburst episode, the low metallicity environment, and the 
numerous clusters, make Mrk 930 galaxy of great interest to 
include in our statistical a nalysis of star cluster populations 
in extreme environments l|Ostlin et al.1 120031 ; lAdamo et al.l 
l2010al . lAdamo et alJlioill) . 

A 3-color composition of Mrk 930, obtained with the 
high- resolution HST Planetary Camera (Figure [T}, shows 
a morphologically perturbed galaxy. Nebular emission sur- 
rounds the burst regions (the green regions are caused by 
strong nebular emission lines falling in the band pass of the 
wide V band filter F606W), confirming the youth of the stel- 
lar population. Knot A is the brightest region in the galaxy. 
It is also the most rich in star clusters, resemblin g the clus- 
ter complexes observed in M82 (|Smith et al.ll2006T l and Haro 



1 Value taken from NED, corresponding to a distance modulus 
of 34.27 mag and a recession velocity of 5485±10 km/s. These 
values are used in the present analysis. 
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11 (|Adamo et alj|2010a[ ). In the Southern part, an extended 
tail is observed. The color of the tail is redder than the main 
galactic body, suggesting an older stellar population. Such 
an evolved population has also been detected in the out- 
skirts of the galaxy (Micheva et al. in prep.). The northern 
part of the galaxy is formed by two close cluster-rich regions 
B and C and a tidal tail. Optical Fabry-Perot interferomet- 
ric imaging shows signatures of a recent merger event in the 
velocity field of the galaxy (Marquart et al. in prep. 2011). 

Here, we present a complete multiwavelength analysis 
of the cluster population in Mrk 930. The data are described 
in Section O The cluster analysis and derived properties are 
discussed in Section [3] In this section we propose possible 
scenarios to explain the origin of the red excess. In Section 
0] we map the mean extinction distribution in the galaxy, 
using ground-based spectroscopy and imaging. We compare 
this mapping with the extinction derived locally from the 
clusters. In the following Section [S] we discuss the proper- 
ties of the galaxy as derived from the cluster formation and 
star formation history (CFH and SFH, respectively). Con- 
clusions are summarized in the final section. 

Throughout the paper we will use the Vega magnitude 
system. 



2 THE DATA 

2.1 HST data: Observations & Photometric 
analysis 

HST multiband high-resolution imaging of MRK 930 were 
carried out in 2007 (associated with program # GO 10902, 
PI: G. Ostlin.). The galaxy was sampled w ith the same 
data set as for the ESO 185-IG13 target (|Adamo et al.l 
|2011| ). The ACS solar blind channel (SBC) camera was used 
to image the galaxy in the FUV filter, F140LP. The op- 
tical sampling of the galaxy was originally scheduled for 
ACS/HRC and WFC but due to the failure of the in- 
strument in 2007 we switched to WFPC2 and in partic- 
ular the Planetary Camera (PC) aperture and obtained 
mages in the U (F336W), B (F439W), R (F606W), and 

1 (F814W) bands . Finally, the IR (H band, F160W) imag- 
ing was performed by the NICMOS3 camera, which offers 
lower spatial resolution but wider field of view. The frames 
were r educed, drizzled and ali g ned using the MULTIDRIZ ZLE 
task (|Koekemoer et all |2002| ; iFruchter fc Hookl 120021 ) in 
PyRAF/STSDAfl The WFPC2 F606W and F814W, the 
SBC/F140LP and the NIC 3 imaging were achieved using 
a dither pattern to improve the final resolution of the sci- 
ence frames. In the case of the F336W and F439W imaging, 
two exposures each were taken to perform cosmic ray cor- 
rection, but not dithering was applied due to the shorter 
available exposure time. The final UV and optical frames 
were rescaled to 0.025" /pixel. The best resolution we could 
achieve for the NIC3 final frame was of 0.067" /pixel. Aper- 
ture photometry on the sources was done using a radius 
of 0.1". The sky annulus around the source was placed at 
0.125" and had a width of 0.05". Due to the crowding of the 
regions, we preferred to use the same aperture radius and 

2 STSDAS and PyRAF are products of the Space Telescope Sci- 
ence Institute, which is operated by AURA for NASA 
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Figure 2. Color magnitude diagram of the cluster population. 
Top panel: the small green (grey) asterisks show the 477 cluster 
candidates detected in the galaxy. The black filled dots are the 
207 clusters with photometric errors in / and R, cr m ^ 0.2. The 
extinction vector indicates the shift of a cluster in this plot, if 
corrected for E(B-V)=0.3 mag. The evolutionary track from Z01 
are for Z=0.004 and several values of the mass. In each model 
track, the squares, diamonds, and triangles indicate 1 Myr, 10 
Myr, and 14 Gyr, respectively. See main text for details. 

sky annulus in all the frames and correct for the missing 
flux (caused by a fixed aperture), using estimated aperture 
correction values for each frame (see Table [TJ . A detailed 
description of the data reduction, cluster catalogue extrac- 
tion, point source ph otometry and charge t r ansfe r efficiency 
correction is given in lAdamo et alj (|2010al . l201ll ) . A list of 
the filters, total exposure time, zero points (ZPs), and other 
observational properties related to the reduction of the data 
are summarized in the Ta ble [T] 

As already found in lAdamo et all (|2010al . 120111 ) the 
number of clusters detected in the two deepest exposures, R 
and I bands, is greater than in the other available filters (see 
Table [1} . In total we recovered 477 cluster candidates (with 
a ^ 1.0 mag). Viewing the color-magnitude diagram (CMD; 
Figure [5J) and evolutionary tracks we expect the cluster pop- 
ulation to show a wide range of masses and ages. After a 
photometric error cut at a 0.2 mag, 207 objects remain. 
Of these 157 sources are detected in at least a third filter, so 
this is the number of objects we consider in the analysis of 
Mrk93C0. The final photometric catalogue contains cluster 
candidates that have been detected in at least three filters 
with a signal-to-noise ratio of 5 or larger (i.e., photometric 
error of a m < 0.2). In Table [T] we show the corresponding 
magnitude limits and the number of objects retained in each 
frame after this selection. The final error associated with the 
integrated flux in each filter is the combined quadrature of 
the photometric error and the uncertainty on the aperture 
corre ction determination . The foreground galactic extinc- 
tion (|Schlegel et al.|[l99ct ) correction is also included. 

The completeness test was performed on the two R and 
I WFPC2 frames, the deepest filters used to extract the 
first cross-correlated source catalogue (|Adamo et al . 2010a. 

3 The complete photometric catalogue is available on request 
from the authors. 
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Figure 1. Three-color image of the starburst galaxy Mrk930: WFPC2 F814W filter in red, F606W in green and F336W in blue. Some 
intense blue spots are caused by an imperfect removal of cosmic rays from the WFPC2/F336W frame. Different star forming regions can 
be observed in the galaxy. The youngest are surrounded by intense nebular emission (green diffuse light, due to the wide V band filter 
F606W which transmits [Oin] and Ha). Knot A is the brightest optical and UV region in the galaxy and is undergoing the most active 
starburst event. Region B and C are less clustered and coeval to knot A. The color of the extended tail in the South of the galaxy is 
redder than the starburst regions, probably formed by a more evolved stellar population. See the main text for more details. 



Table 1. HST observations carried out within the framework of program # GO 10902 (PI G. Ostlin). For each filter we list in the table: 
the total exposure time; the corresponding ZPs (Vega magnitude system); the aperture corrections, a c ; the number of objects detected 
with a S/N ^ 5; and the corresponding magnitude limit, a) Both the HST filter names and the abbreviated nomenclature indicated in 
bracket are used thereafter in the text. 



Instrument 


Filter" 


Camera 


Exposure time 


ZP (mag) 


a c (moj) 


N(cr s; 0.2) 


mag limit 


ACS 


F140LP (FUV) 


SBC 


2532 s 


20.92 


-0.54±0.05 


117 


24.0 


WFPC2 


F336W (£7) 


PC 


1200 s 


19.43 


-0.42±0.05 


81 


23.0 


WFPC2 


F439W (B) 


PC 


800 s 


20.88 


-0.36±0.05 


79 


24.1 


WFPC2 


F606W (R) 


PC 


4000 s 


22.89 


-0.61±0.09 


207 


27.1 


WFPC2 


F814W (J) 


PC 


4500 s 


21.64 


-0.73±0.04 


207 


26.0 


NICMOS 


F160W (H) 


NIC3 


4992s 


21.88 


-2.45±0.37 


78 


25.9 
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Figure 3. Completeness fraction as function of the object 
magnitudes. The outputs of the F606W frame are plotted in 
green(grey), while in black wc show the fraction for the F814W 
frame. See the main text for the definition of "inside" and "out- 
side" regions. The red dotted straight lines show the 90 % and 50 
% completeness limits. 

l201lf ). This type of analysis is useful for constraining the 
sensitivity of our method to recover objects. Similar to the 
Haro 11 and ESO 185-IG13 analysis, we defined a center of 
the galaxy (a = 23 : 31 : 58.5; S = 28 : 56 : 52.0) and a 
crowded region of 4.5" radius, where most of the starburst 
knots were contained. The surrounding region (an annulus 
of inner radius of 4.5" and outer radius of 12.25") is char- 
acterized by a low sky background, ideal for detecting faint 
objects. The number of recovered objects was compared to 
the crowding of the considered region. In Figure [3j we show 
the total fraction of recovered detections. The two solid lines 
(green and black) show that in general we are able to de- 
tect most (90 %) of the objects with luminosity brighter of 
~ 26.8 mag. At roughly 28 mag we are still able to detect 
50 % of the objects in both filters. However this is an aver- 
aged result. Looking at the recovered fraction in the "inner" 
(crowded) and "outer" (starburst outskirts) areas, one can 
see that the detection of objects drops more rapidly in the 
crowded regions and changes quite abruptly. This behaviour 
is caused by the small area of the starburst region which 
does not allowed to build a statistically populated sample 
of objects at different luminosities bins. Since the quality 
and sensitivity of the used data for the three BCGs are sim- 
ilar (same exposure time and observational strategy), we can 
use the results of the completeness tests of Haro 11 and ESO 
185-IG13 to consider approximately 26 mag as a threshold 
for 90% completeness detection in the inner region for both 
filters. 

2.2 Ground based imaging and spectroscopy 

2.2.1 NOT Narrow-band imaging 

Ground based narrow band imaging was carried out with the 
2.56m Nordic Optical Telescope (NOT) on La Palma with 
the ALFOSC instrument. The Ha observations were carried 
out in Sept 2002. We obtained 1800s of integration through 
NOT filter #70 approximately centered on Ha at the red- 
shift of Mrk 930. For continuum subtraction we obtained 



2400s of integration through filter #78. We also obtained 
exposures in both filters of the spectrometric standard stars 
BD+28d4211 and GD248. The data were reduced in a stan- 
dard manner (bias, flatfielding, zeropoint calibration). To 
obtain the Ha flux in the continuum subtracted image we 
applied the following additional steps: i) We corrected for 
the relative transmission of the filter at the redshifted wave- 
length of Ha. ii) We corrected for the contribution of [N n] 
by taking the relative filter transmission at the redshifted 
wavelengths into account and by using spectroscopically de- 
termined line ratios which we assumed to have a spatially 
constant value over the face of the galaxy (this is reasonable 
due to the similar ionization potentials), iii) We multiplied 
with the equivalent rectangular width (the integral of the 
filter throughput curve divided by the filter peak transmis- 
sion) of the filter to obtain an image calibrated in units of 
erg/s/cm 2 . 

The H/3 observations were carried out in August 2007, 
using filters #113 (t oxp = 4900s) to capture H/3 and #17 
(Str6mgren-&, t cxp = 2100s) for the continuum, and the 
spectrophotometric standard FeigellO for calibrations. The 
same reduction and calibration steps as for Ha were followed 
except that we did not need to correct for [Nil]. 

In addition we obtained exposures in filters #66 (600s) 
and #18 (Stromgren-y, 600s) to capture the redshifted 
[Oiii]a5007 line and its continuum. Seeing was close to one 
arcsecond for all observations and conditions were photo- 
metric. 



2.2.2 Spectroscopy 

Low-resolution long slit spectra of Mrk 930 along two po- 
sition angles (PAs) were taken on November 14, 2009 with 
the 3.6m Telescopio Nazionale Galileo (TNG) using the DO- 
LORES focal reducer. We used the LR-B grism which yields 
a wavelength coverage between 3000 and 8430 Aand a resolu- 
tion R~600 for our slit width of 1". The data were acquired 
in 1x2 pixel binning, resulting in a spatial scale of 0'.'504 per 
pixel. Several spectroscopic standard stars were observed at 
various airmasses (AM) for the sake of flux calibration. The 
observations consisted of three consecutive 12 min exposures 
along each position angle and were carried out in either case 
at a low airmass (0.03). No corrections for differential at- 
mospheric refraction were therefore applied. 

Conditions were relatively good during the first point- 
ing with the slit rotated to PA=74.1 deg and encompass- 
ing regions B and C; comparison of multiple science and 
calibration exposures suggest transparency varied by <5%. 
Despite the mediocre seeing (l'.'5 to l'/9 FWHM), we could 
clearly separate regions B and C. For the other slit position 
the seeing was worse, «2'.'5, rendering a spatially resolved 
study of component A (PA=108 deg) impossible. An inter- 
comparison of individual galaxy and standard star spectra 
taken during that phase of the observing run has revealed 
transparency variations by a factor of two, resulting in a sig- 
nificantly lower signal-to-noise (S/N) ratio for at least one of 
the three spectra taken. The data were reduced with MIDAS 
using a standard procedure. 
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EXTRACTING THE CLUSTER 
PROPERTIES FROM THEIR LUMINOSITY 



3.1 The Cluster Luminosity Functions 

It is known that young star clusters have power law lumi- 
nosity (CLF, dN(L) oc L~ a dL) and mass functions (CMF, 
dN(M) = CM' 7 dM), with slopes of ~ 2.0. However, the 
CLF varies from galaxy to galaxy, showing a wide range 



of slopes, 1. 



de Griis et al.l 



sg a ^ 2.4 (lLarsenll2002l; iHunter et al l 
l2003t iMora et alj (20091 iPellerin et all 



2003 



2010 



Whitmore et al.ll2010l ; among many others), or even a bend 
occurring roughly between —8.0 and —10.0 mag (Gieles et 
al. 2006a,b). 

The CLF is composed by many cluster populations, 
formed in different periods of the star formation history of 
the host. Absence or presence of variations in the slope of the 
CLF could probe: i) whether the CMF is a single power law 
with index —2.0 fsee lChandar et "ai1l2010l . IWhitmore et al.l 
I2010I ) or has a Schechter function with a "soft" truncation at 
a characteristic mass, rela ted to the environmental proper- 
ties of the ho st galaxy (see lGieles et al.ll2006bl . lBastianll2008l . 
iLarsenl [20091 ). ii) whethe r cluster disr uption is mass depen- 
dent or independent (see lLamersll2009l for a short summary). 

In studies of star cluster properties in BCGs it has been 
argued that the CLF in these system s is shallower tha n 
expected: a ~ 1.8 in ESO 338-IG04 (|6stlin et alj|2003h ; 
q ~ 1.5 (but a ~ 1.6 if only the brightest b ins are included 
in the fit) in Haro 11 (| Adamo et alj|2010a[ ): a ~ 1.8 (but 
a ~ 2.2 including on ly the brightest bins) in ESO 185-IG13 
(|Adamo et alj|201lh . 

We present here the analysis of the CLF in Mrk 930. We 
use the numerous objects detected in the two deepest expo- 
sures, F606W and F814W, with photometric errors below 
1.0 mag, to determine the slope of the CLF (similar results 
have been recovered if only the a ^ 0.2 sample is used). In 
the left and central panels of Figure f3] we show the two re- 
covered CLFs corrected for completeness. In the right panel 
we include the CLF for the H band. In this last case, it 
is important to mention that we included all the positive 
detections obtained, making the derived distribution rather 
uncertain. 

The fit to the R and / CLFs, showed as a black solid 
lines, includes all the luminosity bins up to the value where 
we recovered 90 % of objects in the completeness tests. The 
fit to the H band CLF excludes the fainter bins where we 
notice a decrement in the number of objects. The slopes 
in the three filters are, inside the uncertainties, similar and 
much flatter than expected. A second fit, confined to the 
brightest luminosity bins, produces higher values, slightly 
more consistent with the expected 2.0 value. 

W e have already discussed in lAdamo et al. | (|2010al . 
l201ll ) that blending can flatten the observed CLF. More- 
over, the incompleteness could be much more severe than 
the values obtained in this analysis. Therefore, it is not pos- 
sible to derive any conclusion with the current data. Simu- 
lations, which reproduce blending effect and different CMF 
shapes, need to be used to explore the physical properties 
of the host environments by means of the CLF. 



3.2 Analysis method: x fitting 

To extract the physical properties of the clusters, i.e., age, 
mass, and extinction, we compared spectral synthesis models 
to their reconstructed spectral energy distributions (SEDs) 
from multiband photometry. The models we used contain 
several assumptions and parameters: 

i) clusters form in an instantaneous burst, i.e., the stellar 
population is coeval; 

ii) the stars in the cluster have m asses sampled by a 
Kroupa universal initial mass function (|Kroupall200lh : 

iii) both photoionized gas (important during the first 10 
Myrs of the cluster evolution) and stars contribute to the 
total integrated fluxes; 

iv) standard H 11 region val ues for the gas as inp ut parame- 
ters to Cloudy (version 90.05. [Ferland et al.lll998T ): hydrogen 
gas density of 10 2 cm~ 3 , filling factor of 0.01, covering fac- 
tor of 1 (there is no-leakage of ionizing photons, e.g. all the 
Lyman continuum photons ionize the gas); 

v) the metallicity of the stars and the gas in the cluster 
are the same, i.e., Z = 0.004 (corres ponding to the values 
determined bv llzotov fc Thuan|[l998l ') . 

A deta iled description of th e mod els can be found in 
Z01 and in I Adamo et all (l2010al. l201ll ). The spectral syn- 
thesis models of Marigo et al. ( 20081 . hereafter M08) have 
been used to test the robustness of some of our results. 

To constrain the model which most closely reproduced 
the observed cluster SED we performed a Y 2 -fit. Applying 
the Calzetti attenuation law (ICalzetti et al. I2000T I. we cre- 
ated a grid of models with values of extinction, E(B — V), 
ranging from 0.0 to 3.0 for each age step. The reddened 
model with the smallest weighted residuals was considered 
the best one and produced age and extinction for the clus- 
ter. The mass was proportional to the normalization fac- 
tor between the best model and the observati on. A detailed 
descrip tion of the x 2 - a lg° r ithm is given in Hd amo et all 
(|2010ah . 



3.3 Reconstructing the SEDs of the clusters 

A close look at the x 2 ~fit outputs reveals that also in 
Mrk 930 there is a f r action of clusters with a NIR excess 
(|Adamo et al.ll2010al . I 201lh . i.e., an observed luminosity at 
wavelengths larger than 0.8 /im which cannot be reproduced 
by our models. The fit procedure for clusters affected by a 
red excess and their analy sis has widely been discussed in 
lAd amo et al. I (|2010al . l201lh . Here, we summarize the main 
characteristics. 

In sources affected by a red excess the fit to the whole 
SED (from UV to IR) fails, e.g. producing high residuals and 
wrong estimates of the cluster age (mass). To trace which 
of the clusters are affected by a red excess, two new fits are 
performed. In these two runs we exclude once the IR band 
(H band) and in the other the I and H bands (the latter 
one if available). The requirement of at least 3 data points 
available, in order to perform the fit, is always retained. 
The outputs of the three fits are then compared and, for 
each cluster, we consider as best fit the one which better 
reproduces the UV and optical shape of the observed cluster 
SED. In general, we notice that even if the I band excess is 
less evident than the IR one and could be easily overlooked, 
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it h as an important impact on t he estimated ages and masses 
(see lAdamo et al.|[2oTOal . l201lT ) . 

In Mrk 930 we performed a \ 2 analysis on a total of 
157 objects. Of these objects, 61 % have normal SEDs, and, 
inside the photometric uncertainties, the performed analysis 
produced fairly good fits. 12 % of clusters display an H band 
flux excess (IR excess) and the remaining 22 % have a flux 
excess at A ^0.8 fim (NIR excess). For these clusters, the 
age, mass, and extinction have been estimated from the UV 
and optical SEDs. In a small number of cluster candidates 
(7 in total, ~5 %), the fit procedure failed therefore they are 
excluded from further analyses. 

Finally, we carried out an individual analysis for a 
rather bright and isolated point-like system, whose nature 
remains unclear. This object is easily noted in Figure Q] on 
the right (west direction) side of Region B, located at a pro- 
jected distance of roughly 1 kpc. It looks reddish, and dif- 
ferent from the color of the clusters in the starburst regions. 
Interestingly, it is not detected in the FUV and U bands 
and is only marginally detected in the B band. Different 
scenarios can be advocated to explain the nature of this ob- 
ject: either a deeply embedded and very youn g cluster (sim- 
ilar c ases have been found in ESO 185-IG13, Adamo et al.l 
l201ll ). or a background object (early type galaxy at redshift 
~ 0.1 or Lyman break system at redshift > 3.0). 

Since it is isolated we performed a new photometric 
analysis, using an aperture radius of 1.0". The model fit 
to the SED of this cluster candidate produced an age of 
35.0 Myr, a mass of 1.2 x 10 7 M o , and E(B - V) = 0.68 
(Av ~ 2.5 and Ab ~ 3.45 mag). The properties constrained 
for this object are quite extreme. The mass is similar to 
the values found for W3 and W30 in the merger remnant 
NGC 7252, and the g lobular cluster G11 4 in NGC 1316 
|Maraston et all 12004 iBastian et alj|2006t ). These evolved 
clusters are located a few kiloparsecs away from their hosts, 
have ages between 400 Myr and 3 Gyr, and dynamical 
masses between 1 — 8 x 1O 7 M0. The formation mechanism 
for such massive systems i s unclear. Numerical simulations 
jFellhauer fc Kroupall200"5h suggest the possibility for these 
systems to form from merger of super star clusters born in- 
side the same cluster complex. 

The cluster-like object in Mrk 930 is much younger than 
these super massive evolved clusters, and even if it poten- 
tially could be a predecessor of such rare objects, we cannot 
exclude any scenario with the current data (i.e., background 



object). The extracted spectrum at the position of this clus- 
ter shows a few faint emission lines (e.g., Ha and H/3) at the 
same rest frame position of the galaxy. The continuum has a 
low S/N and disappeared around 4200 A, in agreement with 
the reconstructed SED from the photometric data. How- 
ever, we note that locally the region surrounding the cluster 
has similar nebular emission lines suggesting contamination 
from the galaxy. Due to the uncertain nature of this cluster 
like object, we exclude it from the cluster analysis. 

3.4 Cluster photometric properties 

Color-color diagrams represent useful tools for analysing the 
numerous Mrk 930 cluster population. In Figure [5] we plot 
3 different colors against R — I, which is used as reference. 
Generally, a negative R — / < color indicates ages younger 
than 10 Myr and is produced by a strong nebular contri- 
bution to the R filter, which includes the bright Ha line. 
However, this assumption is not always valid due to the flux 
excess in the I band. The clusters, in the diagrams, have dif- 
ferent symbols depending of their observed SEDs: black dots 
indicate normal SEDs, red triangles clusters with an excess 
at A ^ 1.0 /im (hereafter, IR excess), and blue diamonds 
SEDs which deviate at A ^ 0.8 /J,m (hereafter, NIR excess). 
We notice that clusters with a NIR excess have a R — I color 
between 0.2 and 1.2 mag redder than the prediction made 
by the best SED fitting model. These residuals can be ob- 
served in Figure |§] and will be discussed in Section f3. 61 The 
R — I color of the "blue" cluster-diamonds is such that, if 
overlooked, causes age (and mass) overestimates, affecting 
the results of the optical based cluster analysis. 

The UV — R (upper left panel) color shows that the 
clusters detected in the FUV are young, at least younger 
than 30-40 Myr. Many of the clusters with an NIR excess 
are located in an area where R — I > 0, e.g., with ages older 
than 10 Myr. The FUV band is sensitive to the reddening. 
Looking at the color-color diagram one can see that clusters 
detected in the FUV have extinctions, E(B - V) < 0.3 ( 
e.g., the arrow in the plot). Therefore, we consider these very 
young FUV bright clusters as systems which have already 
gone through the deeply embedded phase. 

A combination of optical colors (upper right panel) 
probes the impact of the contribution from photoionized gas. 
We compare our models ( solid black line) to the Padova 
tracks (dashed black line: iMarigo et aT]|2008r ). These last 
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Figure 5. Color-color diagrams of the cluster population in Mrk930. Different filter combinations are compared to the R — I color 
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M08 tracks (dashed lines) as well. The black triangles show the starting point of the tracks (1 Myr for Z01 and 4 Myr for M08), the 
red asterisks the position in the tracks at 10 Myr, and the black squares at 14 Gyr. The black dots are clusters with fully fitted SEDs 
(from UV to IR). The red triangles are cluster with an excess at A > 1.0 fj,m. The blue diamonds shows objects with an excess starting 
longword A ^ 0.8 /an. The black arrows indicate in which direction and quantity the colors of the clusters change if a E(B — V) = 0.3 
is applied. Mean errors are included. The green square shows the position of the cluster-like system in the color-color diagram, and the 
green arrow the position after an extinction correction of E(B — V) = 0.68 is applied (see Section I3.3D ■ 



models offer a better modelling of the asymptotic giant 
brach (the difference is evident in the IR color, bottom 
panel). However, due to the age range of the cluster pop- 
ulation we are studying, the impact is negligible compared 
to having self-consistent treatment of the contributing nebu- 
lar luminosity. More than 50 % of the clusters in the optical 
color are located in areas of the diagram that can only be 
reproduced by models that include the photoionized gas. 

Finally, the inclusion of the IR color (bottom panel), 
clearly show that clusters with a flux excess at the redder 
wavelengths are mainly located in an area of the diagram 
where the R — I > and H — R> 1.0 mag. The IR color of 
the clusters with a red excess suggests that the extinction in 
these objects should be much higher than the one predicted 
by the optical and UV colors. In other words, at the NIR 
wavelengths it is possible that we are probing a different, 
more deeply embedded stellar populations, indicating that 
these are very young clusters. We will discuss the cluster 



properties and possible connections with the origin of the 
red excess in Section T3. 61 



3.5 Age, mass, and extinction of the cluster 
population 

In this section we present the derived age, mass, and extinc- 
tion distributions of the cluster population. For 57 % of the 
systems we had at least 4 different data points (integrated 
fluxes) to perform the SED fit. The remaining fits were ob- 
tained from only 3 data points: 13 % with R, I, and H fluxes 
(hereafter, RIH fit); 20 % with FUV, R, and I (hereafter, 
FUV - RI); 10 % with a detection in B or U and R, and /. 

The uncertainties associated with the recovered values 
depend on the used filter combinations; the photometric er- 
rors (including uncertainties produced by the quality of the 
data and the reduction process); the models (limited by our 
ability to reproduce the luminosity properties of the clus- 
ters, by assumptions on the initial mass function - IMF - 
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Figure 6. Final age, mass, and extinction distributions for 150 
clusters in solid black line histograms. The plots include also the 
individual distributions of the 20 clusters detected in R, I, and H 
bands (dotted thick green line) and the 31 clusters for which prop- 
erties were found using FUV, R, and / detections (solid purple 
line). See the main text for details. 



and on the mass-to- light ratio); and the applied attenuation 
law. 

In the cluster analysis, an important factor of uncer- 
tainty is produced by stochastic effects in the sampling of the 
stellar IMF. Se v eral works (e.g. ICervino fc Luridianall2004l : 
iMaiz Apellirhl 120091 ; ISilva- Villa fc LarserJ 1201 ll , among 



many others) have pointed out that stochasticity becomes 
important in low mass clusters (< 10 5 Mq) and may drasti- 
cally affect the photometric studies of clusters with masses 
lower than 1O 4 M0. The spectral evolutionary models of star 
clusters are built assuming a continuously populated IMF. 
In reality, in low mass clusters the stars randomly sample the 
IMF, meaning that they may have formed a number of mas- 
sive stars different from the predicted ones by the models. 
Massive stars dominate th e colors of the clusters. Recently, 
ISilva- Villa fc Larsenl l|201ll ) have showed that the luminosity 
properties of the clusters of < 1O 4 M0 change so drastically, 
that the uncertainties in the age can be as big as 2.5 dex. 
A wrong age estimate has important consequences for the 
mass as well (e.g. if the age is overestimated, the mass will 
also be bigger than the actual value because of the increas- 
ing mass-to- light value as function of the age). However, 
we show in Figure [6] and [7] that many clusters have masses 
higher than 1O 4 M0 suggesting that the stochastic effects, 
although present, are not very important (the error in th e 
age estimate is about 5-10 % ICervino fc Luridianal 120041 ). 
In Section fe.GI we again address stochasticity to explain the 
NIR excess observed in some clusters which luminosity could 
be do minated by red super giant stars (RSGs). 

In lAdamo et alj ()201lh . we discussed the uncertainties 
introduced by the set of data on the recovered estimates 
using Monte Carlo (MC) simulations. Since the data set, 
the models, and the extinction law applied are the same 
as the ones used in the analysis performed on Mrk 930, we 
apply those recovered uncertainties to the present results. 

The recovered age distribution shows, as expected, a 
very young cluster population (upper panel in Figure [B]). 
The most populated age bin is the 3-4 Myr one (log(age) = 
6.5 Myr), similar to what already observe d in Haro 11 
dAdamo et al.l l2010aT ) and ESQ 185-IG13 (|Adamo et all 
l201lf ). From MC tests |Adamo et al.ll201ll ) we know that 
the robustness of this peak is quite high (this age is recov- 
ered in 99 % of the cases). 

A secondary peak is found at slightly older ages, be- 
tween 10 and 20 Myr (signed with vertical black dashed 
lines, log(age)= 7.1 and 7.3 Myr). These two age bins are 
the most uncertain in our analysis. Independently of the 
filter combination used, we observed in the MC tests that 
the recovered number of systems in the 10-20 Myr bins are 
doubled with respect the initial number of objects, mean- 
while the age range between 20 and 150 Myr are emptied. 
The causes of this behaviour is the loop in the evolutionary 
tracks at ages between 10 and 40 Myr (see upper right and 
bottom panels in Figure [S| and a degeneracy between age 
and extinction between 40 and 150 Myr, i.e. the extinction 
vector lies roughly parallel to the evo lutionary tracks (see 
also Section 3.3 in I Adamo et al.ll201ll ). Estimates based on 
a 3-data point SEDs are, in general, even more uncertain. 
The FUV — RI fit produces two spurious peaks at 12 Myr 
(log(age)= 7.1) and at 1 Gyr (log(age)= 9.1). The age dis- 
tribution of clusters with a FUV — RI fit (purple solid line 
histograms, upper panel in Figure[6]) has a considerable frac- 
tion of objects in the 12 Myr bin. Likely, two thirds of these 
objects are in reality older. Similarly, the RIH fit produces 
spurious peaks at log(age)= 6.7 (dash-dot line in Fig. [6| 
and log (age) = 7.3 Myr. We see, however, that the num- 
ber of objects in these two bins is small. We conclude that 
likely the secondary peak is not securely established and 
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Figure 7. Mass-age diagram. Clusters which do not show any 
red excess are represented with filled black dots. Clusters affected 
by IR excess are shown with red triangles. The blue diamonds 
are for clusters affected by NIR excess. The lines (see the insert) 
show the detection boundaries in ages and masses corresponding 
to the magnitude limits reached at a ^ 0.2 in each filter. As 
already expected there are no clusters sitting below these limits. 
The green square shows the mass and age of the cluster-like object 
presented in Section 1331 

open for a scenario where cluster production in the galaxy 
has started roughly 100 Myr ago and continued until the 
present time. Two much older peaks at 1 and 14 Gyr prove 
that the present starburst phase is not the first burst event 
in the galaxy. 

The cluster mass distribution (middle panel in Figure[6]) 
covers a wide range of masses, from 10 3 M© to ~ 10 6 Mq. 
We observe that more than 20 % of clusters are as mas- 
sive or more massive than 1O 5 M0 (the super star cluster 
range). MC tests showed that a significant fraction (roughly 
30 %) of recovered masses are underestimated. The differ- 
ence between input and recovered masses is, however, not 
higher than a factor of 2 (when FUV — RI and RIH fits 
are performed). On ly in the case of a RIH fit we noticed 
jAdamo et al.l 1201 ll ) that in 20 % of cases the mass was 
overestimated up to a factor of 2. 

The mass-age plot (Figure [7| shows the whole cluster 
population. In contrast with Haro 11 and ESO 185-IG13, 
we find no young clusters with masses ~ 10 6 Mq. Clus- 
ters with a NIR excess are all (except one) younger than 
6 Myr. IR excess clusters have ages between 1 and 35 Myr. 
Even in Mrk 930 we observe a lack of low mass (< 10 4 M Q ) 
an d very young cluster s (< 4 Myr). As already discussed 
in lAdamo et al.l (|2010al . l201ll ). a fraction of these objects 
is likely not detected due to blending and crowding. We 
also expect that a considerable number of these clusters 
are still partially embedded in their birth cocoons of dust 
and gas. We estimated that a rather moderate extinction of 
E(B — V) = 0.5 mag (Av ~ 1.5) is enough to make a cluster 
of mass below 10 4 Mq undetected at visual wavelengths. 

The extinction distribution (bottom panel in Figure [6]) 
shows, as expected, that more than 80 % of the clusters have 
extinction E(B — V) ^ 0.2 mag. The remaining fraction is 
distributed in somewhat higher values, reaching 1.2 mag. 
Clusters with extinction higher than 0.4 mag could only be 



Figure 8. Extinction-age diagram. The blue (grey) asterisks show 
each single cluster. Superimposed the mean extinction estimated 
at each model age step and error bars given by the corresponding 
standard deviations at each age bin. 
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Figure 9. Strength of the flux excess in H and / bands as function 
of the cluster properties derived from fitting the blue side of the 
SEDs. The vertical dotted lines in the left side panels separate 
clusters younger than 6 Myr (filled dots) from older clusters (open 
squares). In red we show clusters affected by IR excess (H band), 
and in blue the ones with NIR excess (I and H bands). The 
black crosses show, for each bin, the median value of the residuals 
for normal clusters in the sample, represented by the underlining 
small green ones. 



analysed with a RIH fit. The extinction-age plot (Figure^ 
clearly shows that those clusters with higher extinction are, 
in all cases (except 2), systems younger than 10 Myr. Since 
these clusters are only observed at the redder wavebands 
they are likely still partially embedded. Similar trends in 
the extinction distribution and extinction-age diagram have 
also been observed in the other BCGs, Haro 11 and ESO 
185, and i n systems like the A ntennae (|Mengel et al.l [20051 ) 
and M51 (|Bastian et alj|2005l ). 



© 2011 RAS, MNRAS 000, [111201 



Star clusters in Mrk930 11 



3.6 The red excess in young star clusters 

As already done in the previous studies of Haro 11 and ESO 
185-IG13, we try to trace any potential relation between the 
UV and optical determined cluster properties and the inten- 
sity of the flux excess with respect to the best fitting mod- 
els. In Figure [9] we plot the flux excess (residuals) against 
the corresponding age, mass, and extinction. Represented by 
dots are the clusters younger than 6 Myr and squares older 
systems. The color of the symbols are related to the type of 
excess constrained. If the observed excess is a A > 1.0 lira 
they are in red. If the excess is longword A > 0.8 (im they 
are in blue. As comparison, we include the distribution of 
the "normal" clusters, that are objects with regular SEDs 
(represented by the black and green crosses). 

If we look at the age, mass, and extinction of the clus- 
ters versus their residuals (Am = m mo d — rn bs) in the H 
band (upper panels), we do not see any clear relation be- 
tween the UV-optical properties and the strength of the IR 
excess. We already noticed from the IR-optical color-color 
diagram (bottom panel, Figure [5| that the clusters affected 
by a red excess require much higher extinction values than 
the ones from the UV-optical colors (see the other two dia- 
grams) to move toward the areas of the tracks corresponding 
to the estimated age. This property suggests that another 
mechanism contributes at the IR waveband and not in the 
optical. 

The same analysis, but including the / band residuals 
shows a clear correlation between the age of the cluster and 
the I band excess (bottom panels in Figure [9]). In all except 
one case, the I band is found in clusters younger than 6 
Myr. The spread in mass and extinction is, however, similar 
to the trend observed in the upper panels. 

The age of the clusters (even if determined by UV and 
optical luminosities) may give insights on the origin of the 
red excess. Possibl e origins for t he NIR exce ss have been 
widely discussed in lAdamo et al. l (|2010al . l20Tlb . Here we will 
shortly summarize the main aspects. 

3.6.1 The I band excess 

Due to the narrow cluster age range where the I band excess 
is found, a viable explanation for th is feature is the ext ended 
red emission (ERE, see for a review lWitt fc Viilj|200~j ). The 
ERE is observed as a soft rising continuum between 0.6 — 0.9 
lira. It is observed around galactic and extragalactic Hn re- 
gions (among many others. fperrin et aI.lll99a : lGordon et al.l 
l2000h and caused by a photoluminescence reaction on dust 
grains heated by hard UV radiation. Such energetic photons 
are mainly produced in short-lived massive stars. This could 
explain why the I band excess in our clusters is over after 6 
Myr. A possible contribution from ERE to the / band flux 
of unresolved clusters was also suggested by Reines et al. 
(2008b). 

3.6.2 The IR excess 

Several mechanisms can concur to make the flux at A > 1.0 
(im higher than predicted by models, which include stellar 
continuum and nebular contributions. The distance of the 
galaxy and the resolution achieved - the best with the cur- 
rent accessible facilities - limit our studies to the integrated 



properties of these massive star clusters. However, observa- 
tions of much closer resolved clusters and numerical predic- 
tions of stellar populations in clusters can give us an hint of 
the mechanisms which are likely producing the observed red 
excess. 

Among the youngest and massive resolved star clus- 
ters, 30 Doradus (hereafter 30 Dor) in the Large Magel- 
lanic Cloud (LMC), represents the best reference to un- 
derstand what a recently born star cluster looks like. 30 
Dor is the central region of the extended Tarantula neb- 
ula. Multiwavelength studies of this regions have dissected 
the d ifferent components of th e complex 30 Dor environ- 
ment. IWalborn fc Blades! (|l997h identified five different stel- 
lar populations in the region: the bright core early O-type 
stars which are part of the compact star cluster R136; in 
the north and west region embedded massive YSOs; 3 more 
evolved stellar population groups in the southern and 1.0' 
away in the western region. The R136 cluster has a mas s of 
~ 10 5 A/ o and is 3 Myr old jPortegies Zwart et al]|20ld and 
references therein). This nuclear region (^ 3 pc ) is dust and 
gas- free. Recent studies |Crowther et al.ll2010T ) report that 
stars more massive than the usual assumed theoretical limit 
of 120 Mq reside in the core of R136, probing the extreme 
nature of the star formation process in very massive star 
clusters. However, the aperture radius we are using in our 
analysis is much larger (radius of ~ 36 pc) than the size of 
R136. Our apertures are comparab le to the size of t he im - 
age of 30 Dor showed in Figure 1 of ICampbell et all |2010h . 
This suggests that while the optical range is dominated by 
the stellar and gas emission, the IR transmits also flux from 
the diffuse dust, heated by the hard UV radiation field, the 
embedded YSOs formed in triggered star formation events 
at the edge of the nucleus, where most of the dense gas and 
dust filaments a re located, and the low mass stars still in th e 
PMS phase fsee lWalborn et al.|[2002l ; iBrandner et al.ll200lh . 

In the literature, studied cases of IR excess in young 
embedded or partially embedded unresolved extragalactic 
clusters have explained the red excess as due to an impor tant 
contributi on by YSOs |Fernandez-Ontiveros et all 20091 ') , or 
hot dust (|Cabanac et alj|2005l ; ICresci et all I2005T ). Likely, 
the same mechanism is causing the red excess in young star 
clusters in BCGs. 

After several Myr this complex phase is over, so it can- 
not explain why we still observe objects with an IR excess 
at older ages. For these evolved clusters a possible source 
of excess can be an important contribution from RSGs. 
These are rather massive (~ 25 — 7Mq) stars. They con- 
tribute mainly to the NIR for clusters age between 8 and 60 
Myr. Models usually predict the number of RSGs, assum- 
ing that the stars in a cluster fully populate the underlying 
IMF. However, this assum ption is not valid, if the clu ster is 
less massive than 10 6 Mq |Lancon fc Fouesneau|[2010l ). and 
causes important variations for cluster masses below 10 5 M Q 
(ICervino fc Luridianall2004 ). Moreover, it has been observed 
that in metal-poor environments the number of observed 
RSGs tends to be higher than the predicted one from the 
ratio of blue versus RSGs (see lEggenberger et all [2001 and 
iMaeder fc Mevnetl l200ll ') . Therefore, both effects would be 
observed as a rise in the NIR integrated flux of an unresolved 
cluster, which our current models cannot fully account for. 
It is not clear, however, why we do not see any mass depen- 
dence between the excess in H band and the mass of the 
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Figure 11. The recovered extinction versus the position of the 
cluster in the galaxy. The inset show the color coding used in 
the plot and corresponding to each interval of extinction values. 
The position of the two slits and the regions are indicated. In 
green dots we represented the excluded cluster-like object from 
our analysis. 



cluster, which could support the stochasticity scenario, or 
why our models cannot predict the correct number of RSGs 
only for some of the clusters. NIR spectroscopy is needed to 
prove direct evidences for these scenarios. 

The environment of Mrk930 is quite crowded and we 
cannot exclude that in some cases the red excess could be 
caused by a newly formed cluster in the same aperture or at 
the same line of side, visible only in the IR. 

Other "exotic" explanations for the excess in clusters 
could be a bottom heavy IMF, i.e. a much steeper slope 
at low mass ranges (higher number of low mass stars), or 
a second stellar population formed after a long time-delay 
(~ 10 Myr). Ho wever, studies of IMF variations in mas sive 
star clusters (see lBastian et aj]|20ld : iGreissl et alj|20ich do 
not find any sign of variation; and second populations with 
delay of few tens of Myr have not been observed yet. 



4 MAPPING THE EXTINCTION ACROSS 
THE GALAXY 

We have estimated the extinction in the galaxy using three 
independent analyses: cluster SED fitting, Ha/H/3 line ratio 
from long-slit low-resolution spectroscopy, and narrow band 
Ha/H/3 imaging ratio. 

Regions B and C are resolved in the 2-dimensional spec- 
trum (Pos-1), hence we were able to extract a single spec- 
trum for each of them. The knot A was blended in a single 



blob (Pos-2) so we extracted a single spectrum enclosing 
the total flux. The fluxes of the Ha and H/3 lines were mea- 
sured in the three one-dimensional spectra. To correct the 
two emission lines from the underlying stellar absorption, 
we performed a fit to the continuum using the Starlight 
fitting procedure (see Section I5.3|l . The Ha flux was also 
corrected for the contribution of the blended [N n] line at 
6548.0 A. The flux of this last emission line was estimated 
from the measured flux of the unblended [N n] at 6584.0 A 
and the expected line ratio [Nii16548.o/[Nh16584.o = 0.33 for 
case B recombination (|Osterbrock fc Ferlandl [20061 ). The fi- 
nal extinction was estim ated using the galactic extinction 
law (|Cardelli et al.ll 19891 ). The line fluxes and the estimated 
extinction are all listed in Table [2] 

We constructed an E(B — V) map from the Ha/H/3 im- 
ages assuming a non reddened Ha/H/3 ratio of 2.87 and the 
galactic extinction law. No correction for Balmer absorption 
was applied, and hence the E(B — V) is only indicative and 
may overestimate the extinction in regions where the H/3 
emission is weak. In Figure [10] we show from left to right: 
H/3 image (with logarithmic intensity scaling), and E(B~V) 
map. 

In Figure 1111 we have reconstructed the position of the 
clusters in the galaxy and their recovered extinction. Clus- 
ters with E(B — V) 0.2 (blue triangles and purple di- 
amonds) are the majority. Clusters with higher extinction 
are fewer and located around almost dust-free companions. 
The smaller number of high extinguished clusters could be 
caused by a bias in the selection method. Since optical filters 
have been used to create a position catalogue (see Section 
12. it is likely that those redder clusters have been missed. 
The positions of the two slits are indicated in Figure [TT1 The 
slit is 1" width, however, taking into account the seeing, it is 
likely that the regions close to the edges of the slit have also 
contributed to the transmitted spectrum. In general, we ob- 
serve that the values recovered from the two hydrogen line 
ratios are very low and suggest that the spectra are domi- 
nate by an averaged value of the extinction observed in the 
clusters. 

Similarly, we notice that the resolution of the map in 
Figure [10] is much lower than the one reached with the HST 
data, so the resulting E(B — V) values are averaged over 
a wider region than a typical cluster size. For this reason 
we look at general trends and not to one-by-one analogies 
between a cluster and a pixel-map. Morphologically, we ob- 
serve a correspondence among the brightest H/3 regions (left 
panel, Figure [TO}, the less extinguished areas of the galaxy 
(right panel) , and the location of the very young and, in most 
cases low-extinguished, star clusters (see also Figure I12[) . 
These areas are also the ones which dominate the signal 
of the spectra, because of the very low E(B — V) values 
obtained from the optical line ratios. This feature suggests 
that feedback from young star clusters has likely cleaned 
(destroyed) the dust in these regions (a large scale of the 
30 Dor core). We observe a gradual increase of the extinc- 
tion around the dust-free knots. The centre (between the 
two bright clumps) and the most northern regions of the 
galaxy are the most extinguished. We only detect a dozen of 
clusters in these two regions: many of them have very low ex- 
tinction, or not bigger than 0.4 mag, suggesting that highly 
extinguished clusters in this area have not been detected. 

We conclude that the extinction maps of the galaxy 



© 2011 RAS, MNRAS OOO,rflf20l 



Star clusters in Mrk930 13 




Figure 10. Left: The H/3 image of the galaxy in logaritmic scales (from to 5 X 10~ 16 crg/s/cm 2 ). Ovcrplottcd in white, the contours 
at: 10 — 16 , 3 X 10 — 17 ,10 — 17 ,3 X 10 — 18 erg/s/cm 2 flux levels. The coordinate box units are in arcsec. The color bar is in logarithmic scales 
and goes from f(H/3) value of (black) to 5xl0~ 16 erg/s/cm 2 (white). Right: The extinction map of the galaxy, E(B-V). It is in linear 
scales and the bar cover from (black) to 0.85 (red), the mapped values. The same H/3 contours are overplottod. 



Table 2. Measured f(Ha) and f(H/3) in: the galaxy (first raw) using imaging data, and in the 3 main knots using the extracted spectra. 
The extinction values, E(B — V), are obtained from Ha over H/3 ratio, assuming a Galactic extinction law. The total Ho luminosity, 
L(Hq) has been estimated assuming a distance of 72.0 Mpc. 







f(Ha) 


f(H/3) 


E(B-V) 


L(Hq) 






(10~ 12 erg/s/cm 2 ) 


(10~ 12 erg/s/cm 2 ) 


mag 


(10 41 crg/s) 


Mrk930 


imaging 


1.03 






7.46 


Knot A 


spectr 


0.32 ± 0.08 


0.10 ± 0.17 


0.10 ± 0.04 




Region B 


spectr 


0.23 ± 0.07 


0.07 ± 0.08 


0.16 ± 0.03 




Region C 


spectr 


0.35 ± 0.05 


0.12 ± 0.09 


0.07 ± 0.03 





agree fairly well on the average values, despite having been 
produced with different techniques. Locally, however, only 
the resolution reached by the study of the clusters allow us to 
explore extinction variations on much smaller and detailed 
scales. Of course, there is no a priori reason why the extinc- 
tions derived from nebular lines and continuum should be 
the same even if we could match the HST resolution since 
these spectral components may have different contributions 
from regions along the line of sight. 



5 PROBING THE STARBURST PROPERTIES 
AND THE FORMATION HISTORY IN MRK 
930 

5.1 The bursting nature of Mrk930 

To probe the bursting nature of Mrk930 we constrain the 
^-parameter and the gas consumption time scale, both con- 
sidered starburst tracers. 

The visual magnitude of Mrk930 is My = 15.12 mag 
(Micheva et al. in prep), which corresponds to a total lu- 
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minosity of L v = 3.89 x 10 9 L Q . A more detailed inves- 
tigation of the underlying stellar population will be pre- 
sented by Micheva et al (2011 in prep) but to get an ap- 
proximate estimate of the total stellar mass in Mrk930 we 
take a mass-to-light ratio, M/L= 0.87 which is the me- 
dian for the five luminous BCGs with well constrained stel- 
lar masses in lOstlin et ail (|200lf ). We estimate in this way 
the total stellar mas s in the galaxy, M* ~ 3.4 X 10 9 Mq. 
iHopkins et~aT] (|2002l ) estimated the SFR in Mrk930 using 
the luminosity of the galaxy at 60.0 /im (9.1M0/yr) and at 
1.4 GHz (10.6 M /yr). Estimates of the SFR using the Ha 
luminosi ty give lower values, i.e., SFR(Ha) = 5.34 ± 1.79 
Mq /yr (|Rosa Gonzalez et alj 120091 ) . This value is in very 
good agreement with the estimate we obtained from the 
measured L(Hq) shown in Table [2] Applying the Kennicutt 
relation (|Kennicuttl Il998l ) . we determine a slightly higher 
SFR(Ha)= 5.89 M /yr. 

We constrain an upper limit of the so-called b- 
parameter, i.e. the ratio of the cur rent SFR to the aver- 
age SFR of the past in the galaxy (|Scald Il986l ) . Galaxies 
with constant or declining SFRs have b ^ 1, while systems 
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which are undergoing a violent burst episode have b 3> 1 



which are undergoing a violen 

|Kennicutt et al 1 19941 ). Using the estimates of the M* and 



1400 



of the SFR(Hq), we find b « 16 for Mrk930, confirming that 
the this galaxy is truly a starburst. A similar conclusion can 
be reached from the gas consumption time scale, defined as 
the ratio between the available ga s and the current SF R in 
the galaxy, i.e. T gas =M stls /SFR. iHopkins et al.l (|2002T ) es- 
timated a mass of neutral hydrogen of Mm = 3.02 x 10 9 
Mq . There are not available estimates of the molecular and 
ionized gas mass in the galaxy, so we assume t hat the total 
amou nt of gas is 2 times the determined Mhi (|Lerov et al.l 
120051 ). We find then that Mrk930will be able to sustain the 
current SFR for approximately another half Gyr (r ffas = 1 
Gyr). 

5.2 Starburst propagation revealed by the star 
cluster ages 

To trace the starburst properties and the formation history 
of Mrk930, we compare the cluster properties with other 
probes: the [O m]/H/3 ratio map of the galaxy and the fit to 
the stellar continuum spectra of the 3 starburst knots (next 
sections). 

Using the age of the clusters we mapped the cluster 
formation and propagation in the galaxy (see Figure W2\i . 
We clearly observe that the starburst regions are dominated 
by very young star clusters (blue dots) with an age < 5 
Myr. The [Om]/H/3 ratio is known as a tracer of the mean 
ionisation level and temperature of the photoionized gas 
ijBaldwin et alj|l98lh . The stronger the UV radiation field 
is and the higher the [Om]/H/3 ratio becomes. In Figure [TBI 
one can easily observe that the areas with higher ratio val- 
ues are the ones where the younger star clusters are found. 
Ilzotov fc Thuanl (| 19981 ) estimated for Mrk930 a [Om]/H/3 
ratio of ~ 4.5, in good agreement with an average estimate 
from our map and the classification of this target as a Hn 
galaxy. 

We observe that in the burst regions, very young star 
clusters are located around central older clusters, suggesting 
a progressive propagation of the cluster formation. In the 
extended tail, south of knot A, we mainly find older clusters 
in agreement with the color shown in FigureQ] The tidal tail 
and region C are separated by a group of older clusters all 
located in a small area. The cluster population in the tail is 
rather young and spans an age range of 1-25 Myr. 

Assuming that we are complete in detecting clusters 
more massive than lxl0 4 M o , we derive the CFH during the 
last 40 Myr of star formation in the galaxy. We consider that 
a full population of clusters forms with a CMF power-law 
of index 2 .0 and extrapolate t he missing fraction of cluster 
mass (see lAdamo et alll2010al for details). We obtain then 
a trend of the cluster formation rate in the last 40 Myr (left 
panel in Figure [14)) . which shows a clear increase in the rate 
of cluster formation during the last 20 Myr. The present 
CFR is ~ 1.3M0/yr, i.e., 1.3 Mg/yT of stars are formed 
in bound star clusters. The previous significant burst was 
between 30 and 40 Myr ago. 

Comparing the present CFR to the SFR(Ha), we find 
that ~ 25±10 % of the SF R is happening in bound clusters 
(T =CFR/SFR= 0.25, see iBastianl I2008T ). The fraction is 
quite high, but comparable wit h the values we have f ound 
for Haro 11 and ESO 185-IG13. iGoddard et all |2010T ) pre- 
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Figure 12. Position of the clusters in the galaxy as function of 
their ages (see color coding in the inset). The size of the dots 
is proportional to the estimated mass, i.e., bigger radii are for 
higher masses. The position of the two slits and the regions are 
indicated. The age range between 10-25 Myr (shown in purple) is 
the most uncertain and a fraction of these clusters is likely older 
(< 100 Myr). We also notice that the clusters in the age range 
25-100 Myr (orange dots) are, in reality, all younger than 50 Myr 
(see Figure|6j. In red, we show the few objects detected at roughly 
1 Gyr and in black, the ones at 14 Gyr. See main text for details. 



sented an empirical relation between the cluster formation 
efficiency, T, and the surface density of SFR, Esfr in the 
host. This relation suggests that the SFR and the rate of 
the cluster formation are positively correlated, or in other 
words, that an enhancement in the SFR increases the num- 
ber (mass) of the produced clusters. 

We estimated a radius of the starburst area of the 
galaxy integrating the luminosity flux in B band within in- 
creasing aperture radii, up to include the outskirts of the 
galaxy. We observed that the radius, where 80 % of the light 
produced in the galaxy is collected, included the whole star- 
burst area. We used this radius, r 80 % ~ 1.7 kpc, to estimate 
the Esfr in Mrk930 and found E S fr = 0.55M Q /yr. 

In th e right panel of Figurell4l we show the relation ob- 
tained bv lGoddard et all l|2010l ) and the position of Mrk 930. 
The agreement is good, showing evidence of the a positive 
correlation between the mean properties of the host envi- 
ronment and the production of star clusters. The same cor- 
relation has been found in the others BCGs (|Adamo et al.l 
l2010al . l201lf ). probing that the BCG environment, perturbed 
by the recent merger event, forms a higher f raction of stars 
in bou nd clusters. The results obtained by IWeidner et al.l 
(|2010l ). by means of numerical simulations, support this sce- 
nario. They show that regular patterns, like spiral arms and 
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Figure 14. Left panel: Cluster formation rates during the last 40 Myr of starburst activity. The filled dots connected by the green 
dotted line show the observed CFR derived from clusters more massive than 10 4 Mq. The thick squares show the derived CFR, if the 
total mass in clusters less massive than 1O 4 M0 is extrapolated using a CMF with index —2.0 d own to 10 2 Mp, Right panel: The cluster 
formation efficiency (T) versus the star formation rate surface density, Sgpj^, using data from [Goddard et alj d2010f) (black diamonds 
and the dashed line). The positions of Mrk930 is shown by the green (grey) dot. See text for details. 



Figure 13. The [Oin]/H/3 map of the galaxy produced with 
narrow band imagings. Overplotted in white the H/3contours at: 
10~ 16 , 3x 10" 17 ,10" 17 ,3x 10~ 18 crg/s/cm 2 . The coordinate box 
units are in arcsec. The color bar is linear and goes from [O iii]/H/3 
value of (black) to 10 (red). See main text for details. 



shear, tend to fragment the collapsing GMCs and favour 
the formation of less clustered systems (OB associations). 
On the other hand, in a merging system, the external pres- 



sure on the GMC yields a faster collapse and enhance the 
local star formation efficiency inside the GMC. 

Finally, we notice that in our analyses (Haro 11, ESO 
185-IG13, Mrk930) has not been possible to determine 
which fraction of the clusters we h ave studied are real bound 
cluste rs following the definition bv lGieles fc Porteeies ZwartJ 
1)20101 ). The Goddard et al. relation has been obtained using 
similar analyses. We cannot exclude that an important con- 
tamination by systems which are gravitationally unbound 
could create this apparent relation. Unfortunately, with the 
current data we are not available to test whether the relation 
would still hold or disappea r , hence sugg esting a universal 
cluster formation efficiency (|Bastianll200S ). 



We will discuss this point and the starburst properties 
of the BCGs in a forthcoming paper (Adamo et al. in prep.). 

While the current cluster formation is very high in 
Mrk 930, we find a small number of clusters at age ~ 1 Gyr 
and at ~ 10 Gyr. In the age range around 1 Gyr we find a 
total cluster mass 5.4 x 10 s Mq (clusters more massive than 
10 s Mo). The expected mass in clusters below 10 5 is negligi- 
ble, assuming that clusters in this age range has a Gaussian 
CMF (M peak = 2 x 1O 5 M , g = O.61og(M)M ) . For old 
(age ~ lOGyr) GCs we see a total mass of 4.5 x 1O 6 M 
when including clusters with masses down to 3 x 10 5 Mq. 
Assuming the same Gaussian CMF as in the previous ex- 
ample, the total estimated mass for the old GC population 
is 1.7x 10 7 M o . The stellar mass of the galaxy is M* ~ 3x 10 9 
Mq which would imply that the GCs in Mrk 930 make up 
~ 1 % of the total host stellar mass. T his value is smaller 
than what ha s been found for ES Q 338 (|6stlin et al.ll2003l ) 
and ESO 185 lAdamo et al.l (|201ll ). 
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5.3 Starlight fit to the starburst knot spectra: the 
SFH 

Star clusters are produced in the same global star for- 
mation event which enriches the host of several stel- 



sis codeFl (ICid Fcrnandcs et all 12005 


; ICid Fernandes et al. 


2007 


; iMateus et al.1 20061. Asari et al.1 


2007; see Cairos et al. 



2003 for an applications of Starlight to BCGs) on our spec- 
tra for regions A through C, we try to shed further light on 
the star formation history (SFH) of Mrk 930 and check for 
a consistency the age distribution of stellar clusters and dis- 
crete episodes of enhanced SFR, as revealed from spectral 
synthesis. 

Starlight synthesizes the observed SED of the stellar 
continuum as due to a linear combination of single stellar 
populations (SSPs) of different age and metallicity. The SSP 
library used is based on stellar models bv lBruzual fc Chariot! 
(120031 ) for a Salpeter IMF, three stellar metallicities (Z Q /19, 
Zq/4.75 and Z Q /2A, assuming Z Q =0.019) for 59 ages (be- 
tween 0.25 Myr and 13 Gyr). 

From the best-fitting SSP population vector we derived 
a number of secondary quantities, of which we present here 
the luminosity-weighted and mass-weighted stellar age (see 
Figure [T5"1) . Prior to spectral fitting, several emission lines 
and a spectral region around 7500 A affected by sky line 
subtraction residuals were interactively flagged. 

We remark that the input spectra are not corrected for 
the contribution of nebular continuum emission, thus intro- 
ducing some uncertainty on the outcome of the fitting pro- 
cedure. Ideally, in galaxies with very strong nebular emis- 
sion, the preferable procedure would be to compute and sub- 
tract from the observed spectrum a synthetic nebular line 
and continuum spectrum, in order to recover and separately 
model the spect ral energy distribution (SED) of the stellar 
background (see lPapaderos et al.ll 199^ for an application to 
the metal-poor BCG SBS 0335-052E). However, in the case 
of Mrk 930, since t he EW(H/3) is relatively low (~ 90 A, 
iGuseva et al.|[2000l ). and due to the moderate S/N of our 
spectroscopic data, especially for region A, we prefer not to 
adopt that approach and only flag emission lines. 

As spectral synthesis involves several free parameters 
and the SFH quoted from it is subject to uncertainties, we 
restrict ourselves to a qualitative comparison between the 
main SFH features, as derived from Starlight models, with 
the reconstructed CFH based on our star cluster analysis. 

In the top panels of Figure 1151 we show an example 
of the fit produced by Starlight (the fits to the others two 
spectra are in Figure [T] of the Appendix). The fit to the 
observed spectrum of region B (orange solid line) is fairly 
good, as showed by the residuals in the panel below. An 
hint of the WR-bump observed by llzotov fc Thuanl (|l998T ) 
is also observed in the residual emission spectrum above the 
best-fitting stellar SED. A similar good fit has also been 
obtained for region C, while the outputs residuals of the 
Region A are quite noisy suggesting a slightly poorer fit (see 
Figure [T] of the Appendix). The three plots in the bottom of 
Figure [TS] show the output SFH of the regions as function 
of the fraction of luminosity (top insets) and stellar mass 
(bottom insets). 

4 STARLIGHT & SEAGal: |http://www.starlight.ufsc.br/| 



In region B and C, between 76% and 80% of the ob- 
served flux at the normalization wavelength (5200 A) is 
due to a young stellar population formed within the past 
t y = 10 s yr. The mass fraction of these young stars M y ^ is 
inferred to be 5.5% in region B and up to ~16% in region 
C. 

The dominant luminosity contribution of young stars is 
as well reflected upon the low luminosity-weighted stellar age 
of tL~2.1±l.S Gyr and 0.6±0.3 Gyr for regions B and C, re- 
spectively. This is by a factor of at least 5 lower than the re- 
spective mass-weighted stellar ages tu of 11.5±0.9 Gyr and 
6. Oil. 9 Gyr, suggesting, consistently, that old stars domi- 
nate by mass, more specifically, that 50% of the total stellar 
mass was already at place 10 Gyr ago. With regard to the 
recent SFH of B&C, the bottom plots in Figure [T5] reveal 
at least two salient features: a very recent (<5 Myr) and in- 
tense burst of star formation (SF) which led to an estimated 
4% and 13% of the stellar mass in regions B and C and a 
previous episode of strongly enhanced SF between 10 and 
~40 Myr ago. Spectral fits also suggest that region B has 
been relatively quiescent in the age interval between ~50 
Myr and 1 Gyr, with merely a hint for a minor SF event 
~100 Myr ago. This is not the case, however, for region C 
which, according to our spectral models has undergone sig- 
nificant SF activity between ~200 to 350 Myr ago, with the 
stellar mass formed during that period being comparable to 

My,,. 

The lower S/N ratio of our spectra for region A does 
not permit a conclusive investigation of possible WR fea- 
tures in it. Additionally, due to the poor seeing during the 
observations of that region, it has not been possible to dis- 
entangle and separately fit the spectrum of its NW and SE 
part. Spectral synthesis to the integral spectrum of region 
A indicates that this region experienced a strong burst of 
star formation between 25 and 30 Myr ago which produced 
~11% of its stellar mass, followed by a weaker (~1%) SF 
episode about 5 Myr ago. 



5.4 Comparing star formation and cluster 
formation history 

Spectral fits consistently imply that a strong burst of star 
formation was initiated ~40 Myr ago almost coevally in re- 
gions B and C. This is in good agreement with the analysis 
of star clusters which reveals a sizeable population in the age 
interval between 10 and 40 Myr in either components. Simi- 
larly, spectral synthesis models corroborate the evidence for 
a substantial population of very young (<5 Myr) stars that 
was already drawn from the analysis. 

IGuseva et all i|200d ) estimated from equivalent width 
EW(H/3) measurements a burst age of 4-5 Myr for Mrk 930. 
From the a nalysis of archiv a l Spi tzer Infrared Spectro- 
graph data, lO'Halloran et all |2006T l found for Mrk 930 a 
[Ne ill] /[Nell] ratio of 3.55. The ratio between these two in- 
frared lines is considered an indicator of the hardness of the 
UV radiation field (dominated by the short-lived massive 
stars), and can be converted into an estimate of the burst 
age. Using the estimates produced bv lMadden et all (|2006T ) 
of the [Ne ill] /[Nell] variation as function of the burst age 
(Figure 13 in their paper), we find a consistent upper limit 
for the burst age of 4 Myr. One should be aware, however, 
that the line ratios depend on the ionisation parameter, i.e., 
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Figure 15. In the top panels of Figure 15, we show the Starlight fit for region B (the fits to the other two regions are in Figure 1 of the 
Appendix). The fit (orange solid line) is fairly good, as shown by the residuals in the lower panel. A hint to the WR-bump observed by 
Izotov & Thuan (1998) is also visible in the residual emission spectrum above the best-fitting stellar SED. A comparably good fit has 
also been obtained for region C, while for region A, due to its relatively noisy spectrum, the SFH is more uncertain. The three plots in 
the bottom of Figure 15 show the output SFH of the three regions (A in the left, B center, C right), i.e., age distribution of the SSPs 
selected by Starlight as a function of, respectively, their luminosity contribution (top insets) to the normalization wavelength and mass 
contribution (bottom insets). Thin grey vertical lines in the fraction of luminosity inset indicate the ages of the library SSPs. The shaded 
area in the lower diagram shows a smoothed version of the mass fraction represented by individual SSPs and is meant as schematic 
illustration of the SFH. 



the ratio of the mean photon flux to the mean atom den- 
sity. Thus both the radiation field and the pressure of the 
ambient ISM should be considered. Dynamical modelling 
of the evolution of Hn regions around young clusters show 
that there is a loose coupling between the age and mass of 



the cluster and the ionisation paramete r which then can b e 
constrained to a rather narrow range (|Dopita et al.l 1200(1 ). 
Using the results from Dopita et el. we find that the age de- 
terminations based on the [Nem]/[Nell] line ratios should be 
accurate to within a factor of 2. The age recovered by this 
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line ratio corresponds, however, to a peak in the cluster pro- 
duction, as shown in the previous sections. The agreement 
among different tracers seams to suggest that the burst has 
reached a peak roughly 4 Myr ago, or more in general, that 
it is actively producing stars at the present time. 

There is a considerable amount of gas which can fuel the 
present burst. However, we are not able to predict whether 
this burst episode will evolve into a more quiescent phase or 
will continue converting the gas into stars with a constant 
rate for the next 1 Gyr. 

The Starlight outputs are consistent with the stellar 
component in Mrk930 being predominantly old, with half 
of the stellar mass having being assembled ~10 Gyr ago. 
Despite exhibiting a dramatic burst of star formation at the 
present epoch, Mrk930 seems to be an ancient system. 

We have not found any very old (10-14 Gyr) cluster 
in the starburst regions. On the other hand, we do detect 
some globular clusters with such age range in the galaxy, as 
showed in the mass-age diagram in Figure [7] Their masses, 
however, are a few per mille of the estimated total stellar 
mass of the galaxy. These old clusters are in the outskirts of 
the starburst regions (Figure [T2| where an old stellar popu- 
lation has also been revealed (Micheva et al. in prep). Two 
mechanisms may explain the presence of evolved globular 
cluster in the outskirts and not in the starburst regions. 
Firstly, there are old clusters in these regions but they are 
not detected due to crowding and/or the bright and diffuse 
nebular emissions. Secondly, the merger event has probably 
perturbed the dynamical equilibrium of the involved sys- 
tems. The old galactic stellar population has homogeneously 
settled on the new potential, while the old clusters have 
moved on more distant orbits, similarly to the old globular 
clusters observed at high galactic latitude in the Milky Way. 



starburst region, additional support for the ancient origin of 
Mrk 930. 

The recovered extinction in the very young clusters 
shows considerable spread and decreases at older ages. Our 
analysis suggest that we are mainly seeing the optically 
bright objects, i.e., systems that are still only partially em- 
bedded in their natal cocoons. The deeply embedded clus- 
ters are most likely too faint to be detected at visible wave- 
lengths. We map the extinction across the galaxy using low- 
resolution spectra and ground based imaging (Ha/H/3 ra- 
tio) and we compare these results with the extinction dis- 
tribution derived from the clusters. We find that the mean 
optical extinction derived in the starburst regions is very 
similar to the values typically observed in the clusters (more 
than 80 % of systems have E(B — V) < 0.2 mag), but do 
not trace clusters locally more extinguished. 

Similarly to other blue compact galaxies like ESO 185- 
IG13 and Haroll, a considerable fraction of clusters in 
Mrk 930 shows a flux excess at A ^ 0.8 fj,m, evident from 
our I and H band photometry. The origin of the excess re- 
mains unknown, but we consider a range of possibilities. We 
notice a correlation between the I band excess and the age 
of the clusters (< 6 Myr). The extended red emission re- 
mains the most plausible mechanism which may cause the 
/-band excess. The longer wavelength IR excess that shows 
up in the H-h&nA requires a different mechanism. At young 
ages it could be produced by hot dust and/or a high fraction 
of Young Stellar Objects (hot dust in circumstellar disks). 
Both of these components have been observed in nearby re- 
solved young star clusters. At greater ages, a contribution 
from Red Super Giant stars or other more "exotic" mecha- 
nisms is possible. 



6 SUMMARY 

We have described the starburst properties of the blue com- 
pact galaxy Mrk 930 by undertaking a detailed multiwave- 
length study of its hundreds of resolved star clusters, using 
imaging from the WFPC2 on board the Hubble Space Tele- 
scope, and supported by ground-based narrow-band imag- 
ing. 

Optical and IR line ratios ([Om]/H/3, [Ne ill] /[Nell]) 
suggest that the burst age in Mrk 930 is ~ 4 Myr. The age 
distribution of the star clusters supports this scenario, with 
a prominent peak in the cluster age distribution between 3 
and 4 Myr. The derived cluster formation history shows a 
rapid rise in the last 10-20 Myr, with a rate at the present 
time of 1.3 M@/yr of stars forming in bound star clusters. 
In total, we find that roughly 25 % of the star formation is 
happening in clusters, supporting a scenario where the host 
environment favours the formation of more numerous and 
massive clusters. 

The recent burst episodes traced by the analysis of the 
underlying stellar population are in fairly good agreement 
with the star formation history reconstructed from the clus- 
ters. Young stellar populations dominate the total luminos- 
ity of the galaxy. However, the galaxy assembled more than 
50 % of its total stellar mass more than 10 Gyr ago. We de- 
tect a few, very old globular clusters in the outskirts of the 
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Figure 1. Continue from Figure [TBI Starlight fit to the observed spectra (orange solid line), normalized at 5200 A, of knot A (upper 
panel) and region C (lover one). Superimposed is the best- fitting stellar SED (blue solid line). The corresponding residuals are displayed 
below each panel. The shadowed regions of the spectrum have been excluded from the fit. 
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